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Abstract 

This  paper  presents  results  of  a  joint  research 
effort  of  NASA  Dryden  Flight  Research  Center 
and  the  Institute  of  Flight  Mechanics  and  Flight 
Control  of  Technische  Universitat  MUnchen  on 
longitudinal  flight  dynamics  at  hypersonic  speeds. 

At  high  speeds  a  large  time  lag  between 
pitch  attitude  and  flight  path  angle  exists.  The 
mechanisms  which  cause  this  path-attitude  de¬ 
coupling  are  described.  It  will  be  shown  further 
that  an  limiting  effect  exists. 

An  “aggressive  stick”  piloting  technique  to 
avoid  the  shortcomings  of  the  path-attitude  lag 
is  described.  This  type  of  control  is  also  effec¬ 
tive  to  stabilize  phugoid. 

The  simulation  facility  of  Dryden  Flight 
Research  Center  which  was  used  in  the  evalua¬ 
tion  is  described.  Simulation  results  concerning 
the  effect  of  path-attitude  lag  on  “aggressive 
stick”  technique,  of  time  delay  and  scaling  of 
vertical  speed  information  are  presented.  Fur¬ 
thermore  results  of  flight  tests  with  SR-71  and 
Tu-144  at  high  supersonic  speeds  for  similar 
topics  are  presented. 

Nomenclature 

A(s)  coefficient  matrix  of  homogeneous  system 
Ay  factor  of  flight  path  transfer  function 
Ati  factor  of  pitch  attitude  transfer  function 
B  scaling  matrix  of  control  inputs 


CD  drag  coefficient 
=dCD/d(V/V0) 

CDa  =dCD/da 

CDs  =dCD/d8e 
CL  lift  coefficient 
CLV  =dCL/d(V/V0) 

CL(x  =  dCL  /  doc 
CLS  =dCL/dde 
Cm  pitch  moment  coefficient 

cm  =(1  lp„)dcmldh 
cmq  =2  3Cmld(qclV0) 

C.w  =3  CJd(VIV„) 

c„„  =  3c„,/3a 

c, _  =  23Cm/3<«  c/Vo) 

C.  =  3c„/3r 

my  m  • 

C  _  =  dC  / dS 

mg  m  e 

c  mean  aerodynamic  chord 
g  acceleration  due  to  gravity 
h  altitude 
i  radius  of  gyration 

KEAS  knots  equivalent  airspeed 
M  Mach  number 
m  mass 

Nys  numerator  of  flight  path  transfer  function 
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Ng  numerator  of  pitch  attitude  transfer  function 
nh  thrust  /  altitude  dependence, 

nh=(\/ph)d(T/T0)/dh 
nv  thrust  /  speed  dependence, 

nv  =(V0/T0)dT/dV 

5  reference  area 

s  Laplace  operator 

T  thrust 

T2  time  to  double  amplitude 

Th  height  mode  time  constant 

Tgi  low  frequency  pitch  attitude  zero 

T02  high  frequency  pitch  attitude  zero 

t  time 

u  control  vector 

V  speed 

x  state  vector 

zT  vertical  thrust  line  offset 

a  angle  of  attack 

/  flight  path  angle 

A  denoting  a  perturbation,  e.g.:  AV 

A(s)  denominator  of  transfer  functions 

8e  pitch  control 

8t  throttle  position 

C,p  phugoid  damping  ratio 

C,SP  short  period  damping  ratio 

pitch  attitude  damping  ratio 
p  relative  mass  parameter,  p  =  2  m/{p()  S  c ) 

6  pitch  attitude 

p  air  density 

ph  density  gradient,  ph  -  (1/  p0  )dp  /d h 
<7  damping 

T  reference  time,  z  -  pc/V0 

co  frequency 

(Op  undamped  natural  phugoid  frequency 
c osp  undamped  natural  short  period  frequency 
(O0  undamped  natural  pitch  attitude  frequency 


1  Introduction 

For  improving  performance  and  reducing  launch 
costs  of  space  transportation  systems  a  signifi¬ 
cant  interest  in  aerospace  planes  which  are  using 
aerodynamic  lift  and  air  breathing  engines  exists 
[1-3].  These  vehicles  which  are  capable  of  sus¬ 
tained  hypersonic  flight  require  research  for 
solving  many  challenging  problems. 

Besides  the  fields  of  aerothermodynamics, 
propulsion  and  materials  also  flight  dynamics 
and  associated  handling  qualities  show  unique 
characteristics  [4],  This  is  the  result  of  the  inte¬ 
grated  airframe/propulsion  system  which  is  char¬ 
acteristic  for  aerospace  planes  and  causes  a  sig¬ 
nificant  coupling  between  aerodynamics,  propul¬ 
sion  and  flight  dynamics.  Further  the  high  kinetic 
energy  level  during  hypersonic  flight  and  altitude 
dependencies  have  a  substantial  effect  [5]. 

Right  path  control  is  among  the  unique  prob¬ 
lems  of  hypersonic  flight  dynamics.  In  subsonic 
flight  the  delay  between  pitch  attitude  and  flight  path 
after  a  stick  input  is  small,  usually  below  2  seconds. 
Thus  the  pilot  can  substitute  one  for  the  other.  In  the 
high  speed  regime  a  large  time  lag  of  20  seconds  or 
more  [6]  exists  between  the  pitch  attitude  change 
following  an  elevator  input  and  the  associated  flight 
path  response.  This  reaction  will  be  referred  to  as 
path-attitude  decoupling  or  inconsonance. 

Different  aspects  of  the  path-attitude  decou¬ 
pling  problem  are  subject  of  theoretical  and  ex¬ 
perimental  research  [6-11],  The  difficulty  is  con¬ 
sidered  to  be  associated  with  the  high  speed  [6] 
and  low  lift  curve  slope  [7]  typical  for  vehicles 
capable  of  hypersonic  flight.  Another  reason  is  a 
low  frequency  washout  characteristic  due  to  al¬ 
titude  effects  [8].  Further  flight  controllers  using 
direct  lift  control  have  been  evaluated  [9]. 

2  Dynamics  of  Hypersonic  Flight 

2.1  Equations  of  Motion 

At  hypersonic  speeds  altitude  dependant  forces 
and  moments  are  significant  [5]  and  must  be 
taken  into  account  just  as  engine  effects  [12], 

The  equations  of  motion  in  the  hypersonic 
speed  regime  may  be  described  as 

A(s)  Jc(s)  =B  u(s)  (1) 

where 
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ST 

+  (2  -nv)CD  +  CL v 

CDa 

stT8° 

Vo 

+  li  c  Ph  CD  ( 1  ~nh) 

2  C,  +  CLv 

^La  +  Q> 

-(sT)2+pc  ph  CL 

~(sT)2(iy/c)2 

~(s t)3  ( iy/c )2 

^cmv 

+  ST(Cmq+CmA  ) 

+  (sr)2Cmq  +sTliCmy 

+  ^C>na 

+  V2cphCmh 

(2) 


Jt(j)  =[AV/V0,  a  a,  A  h/(pc)]T  (3) 


CD§ 

-CD 

B  =  - 

CLS 

0 

(  Z>T  !  C  ) 

«<»  =  [§,,  8t]t  (5) 

The  short  term  dynamics  show  similar  de¬ 
pendencies  as  in  the  subsonic  speed  regime.  The 
long  term  dynamics  however  have  unique  char¬ 
acteristics  which  manifests  in  phugoid  and 
height  mode.  The  following  approximations  apply: 


•  Phugoid: 

to2p~-g  ph 

(6) 

(7) 

•  Height  mode: 

\^r  — 

(8) 

'  h 

r  v 

W  V0 

2.2  Path-Attitude  Decoupling 

As  shown  in  Fig.  1  the  time  delay  between  pitch 
attitude  and  flight  path  following  an  elevator 
step  input  in  the  subsonic  speed  regime  is  rather 
small  (c.  Is).  At  hypersonic  speeds  however  a 
large  time  lag  of  20  seconds  or  more  exists  be¬ 
tween  pitch  attitude  and  flight  path  response 


M=0.4;  h=3,000  m 


M  =  1 0 ;  h=34,000  m 


Fig.  1  Time  histories  for  a  pitch  step  input  (-1°  at  1  s)  for  subsonic  (M=0.4,  h= 3,000  m)  and  hyper¬ 
sonic  (GHAME  model,  M=  10,  /?=34,000  m)  flight 
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which  is  referred  to  as  path-attitude  decoupling. 
The  changes  of  the  angles  concerned  are  much 
smaller  in  hypersonic  flight  than  in  the  conven¬ 
tional  speed  regime.  Despite  this  the  resulting 
vertical  speed  however  is  of  the  same  order  of 
magnitude  in  both  cases. 

The  path-attitude  decoupling  effect  also  ma¬ 
nifests  in  the  corresponding  transfer  functions 
which  for  the  frequency  range  of  interest  may 
be  expressed  as 

7(s)_=AjG)=^V 

Se(s)  A (s)  A (s) 

0(s )  _  ATgQ)  _  Ae  Q  +  I/Tqi)  (s  +  1/T92) 

Se(s)~  A(J)  "  A(s) 

where 

A(s)  =  (s2 +2£sp  (osp  s  +  (02sp){s2 +2t;p  (op  s  +  co2 ) 


The  numerator  zero  T$i  shows  different 
dependencies  compared  with  the  subsonic  case. 
It  may  be  expressed  as: 


y0ph 


(id 


Nes(s)  =  Ae(s2  +2£g  cog  s  +  co2g)  (13) 

where 

(°e~-8Ph  (14) 


s  Pa 

2V0  CL\  p 


There  is  no  aerodynamic  or  speed  effect  on  (Oq 
and  it  is  also  similar  to  the  phugoid  frequency  (6). 
At  higher  velocities  the  damping  ratio  converges 
towards  zero.  Therefore  at  higher  Mach  numbers  a 
pole  zero  compensation  between  the  complex  pitch 
attitude  zeros  and  the  phugoid  exists. 

The  described  characteristics  of  the  pitch  at¬ 
titude  zeros  results  in  an  upper  limit  for  the  delay 
between  path  and  attitude  changes  at  a  high  level. 

Fig.  3  depicts  frequency  responses  for  pitch 
attitude  and  flight  path.  Pitch  attitude  has  a  differ¬ 
ent  shape  than  in  the  conventional  speed  regime 
and  is  flat  at  frequencies  between  coP  and  (QSp, 
corresponding  with  path-attitude  decoupling. 
Flight  path  angle  has  a  similar  shape  as  at  sub¬ 
sonic  speeds  but  a  lower  amplitude.  The  region 
between  (Op  and  (Osp  shows  a  KJs  characteristic. 


The  other  zero  Te2  however  is  similar  to 
the  subsonic  case: 


‘02 


8  CLa 


(12) 


As  shown  in  Fig.  2  both  zeros  move  with 
increasing  airspeed  towards  each  other  and 
combine  to  complex  zeros  yielding 


o  [1/s]  -0.04  -0.03  -0.02  -0.01  0 


Fig.  2  Effect  of  Mach  number  on  attitude  nu¬ 
merator  zeros  (GHAME  model,  maxi¬ 
mum  lift-to-drag  ratio) 


j(Op 

j  K  Kegion 

/ 

(SP 

w 

-40 


Fig.  3  Asymptotic  Bode  plots  for  pitch  atti¬ 
tude  and  flight  path  angle  transfer 
functions  (GHAME  model,  M=10, 
h=\  10,000  ft) 
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2.3  Flight  Path  Control  Technique 

For  manual  control  tasks  an  open-loop  K/s  char¬ 
acteristic  in  the  frequency  range  of  concern  is 
preferable  [13].  As  shown  in  Fig.  3  the  flight 
path  angle  response  has  such  a  shape.  Therefore 
it  may  be  used  as  reference  for  flight  path  con¬ 
trol  by  the  pilot. 

Since 

h-V  sin  y  (16) 

vertical  speed  control  is  an  equivalent  reference  for 
flight  path  control  with  speed  approximately  con¬ 
stant.  This  technique  was  investigated  in  simulator 
experiments.  It  was  described  as  “aggressive  stick” 
by  the  involved  pilots.  A  similar  vertical  speed 
feedback  piloting  technique  was  also  used  for  fly¬ 
ing  the  SR-71  [11].  The  piloting  technique  in  mind 
is  described  in  more  detail  in  [14]. 

Fig.  4  presents  a  comparison  between  the 
“aggressive  stick”  and  a  conventional  pitch  con¬ 
trol  technique.  The  vertical  speed  deviations  with 
the  “aggressive  stick”  control  technique  are  much 
smaller  than  those  of  the  other  case.  There  is  also 
a  significant  difference  in  the  Cooper-Harper 
rating  in  favor  of  the  “aggressive  stick”  technique 
(4  vs.  7  in  the  other  case). 

Fig.  5  shows  flight  test  results  for  a  hori¬ 
zontal  turn  maneuver  during  supersonic  cruise 
with  pitch  attitude  (Tu-144)  and  vertical  speed 

.  Conventional  Pitch  Control 


Fig.  4  Comparison  of  conventional  and  ag¬ 
gressive  stick  piloting  techniques 
(GHAME  model,  M=10,  h=l  10,000  ft) 


.  Tu-144  (Pitch  Attitude  Cue) 

-  SR-71  (Vertical  Speed  Cue) 


Fig.  5  Comparison  of  different  pilot  cues  for 
horizontal  turn  task  maneuvers  during 
supersonic  cruise 

(SR-71)  cues  for  the  pilot.  Altitude  and  pitch 
attitude  deviations  with  vertical  speed  control 
technique  are  much  smaller  than  those  with 
conventional  pitch  attitude  control.  This  is 
similar  to  the  hypersonic  results  (Fig.  4). 

An  additional  advantage  of  the  above  flight 
path  control  technique  is  a  stabilization  of  the 
marginally  stable  phugoid.  This  is  shown  in  Fig.  6 
which  presents  the  effect  of  the  flight  path  angle 
feedback  on  phugoid  and  height  mode. 


Fig.  6  Effect  of  flight  path  feedback  to  elevator 
(GHAME  model,  Af=10,  /?=  11  (),()()()  ft) 
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3  Simulator  Experiments 

3.1  Simulation  Facility 

For  the  evaluation  a  fixed  base  simulator  of  the 
Dryden  Flight  Research  Center  with  con¬ 
ventional  stick,  rudder  and  throttle  controls  was 
used  (Fig.  7).  Stick  and  rudder  pedals  are  trima- 
ble  and  the  associated  forces  are  electrically 
generated.  A  Space  Shuttle  type  instrumentation 
panel  is  installed  in  the  simulator.  For  the  ex¬ 
periments  a  customized  simulated  head  up  dis¬ 
play  was  used. 

Test  pilots  of  NASA  Dryden  Flight  Re¬ 
search  Center  with  a  large  experience  in  high 
speed  flight  conducted  the  evaluation. 

3.2  Simulation  Model 

For  simulating  a  hypersonic  vehicle  the 
GHAME  program  [15]  was  used.  It  comprises 
aerothermodynamic  and  engine  models  of  a  ge¬ 
neric  single  stage  aerospace  plane  from  start  to 
orbit  (0  •  M  •  25).  The  six  degree  of  freedom 
equations  of  motion  include  effects  of  an  oblate 
and  rotating  earth. 

Level  I  short  period  handling  qualities  [16] 
have  been  realized  with  an  appropriate  pitch 
damper. 

3.3  Simulated  Head  Up  Display 

The  simulated  head  up  display  (Fig.  8)  was 
customized  for  hypersonic  cruise  tasks. 

In  the  center  of  the  display  is  a  high  reso¬ 
lution  pitch  ladder  moving  around  a  fixed  sym¬ 


bol  (=w=).  Additional  information  on  the 
pitch  ladder  concerns  a  high  resolution  vertical 
speed  indicator  (diamond  symbol  ❖).  Further¬ 
more,  path  and  sideslip  angle  are  presented  (air¬ 
plane  symbol).  Other  scales  indicate  longitudi¬ 
nal  and  vertical  acceleration  as  well  as  speed, 
altitude,  bank  angle  and  heading.  A  more  de¬ 
tailed  description  is  included  in  reference  [17]. 

3.4  Atmospheric  Disturbances 

During  Space  Shuttle  flights  density  changes  at 
constant  altitude  have  been  observed  [18].  In 
accordance  to  those  observations  density  pertur¬ 
bations  of  (1-cos)  type  have  been  implemented 
with  a  maximum  amplitude  of  2.5  %  of  the  ac¬ 
tual  density  and  a  duration  of  20  seconds.  The 
sign  and  interval  between  the  perturbations  is 
selected  by  random. 

3.5  Maneuver 

The  following  representative  maneuvers  for 
hypersonic  cruise  tasks  have  been  specified  with 
corresponding  performance  limits: 

•  Horizontal  turn:  The  task  is  to  perform  a 
heading  change  of  12°  with  constant  equi¬ 
valent  airspeed  and  altitude  followed  by  a 
flight  at  the  final  heading  with  constant 
speed  and  altitude  for  1  minute  (Table  1). 

•  Vertical  plane  altitude  change  :  The  task  is 
to  perform  an  altitude  change  of  3,000  ft 
with  2,000  ft/min  vertical  speed,  followed 
by  a  flight  at  constant  altitude  and  speed 
for  1  minute  (Table  2). 


Fig.  7  Cockpit  of  the  civil  transport  simulator 
used  during  evaluations 


Fig.  8  Simulated  head  up  display  used  during 
evaluations 
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Controlled 

Adequate 

Desired 

Parameters 

Performance 

Performance 

Target  Heading 

+1° 

±0.5° 

Target  Altitude 

±600  ft 

±300  ft 

Trim  Speed  (EAS) 

±10  kt 

±5  kt 

Table  1  Controlled  parameters  and  performance 
values  for  horizontal  turn 


Controlled 

Parameters 

Adequate 

Performance 

Desired 

Performance 

Target  Altitude 

±600  ft 

±300  ft 

Trim  Speed  (EAS) 

±10  kt 

±5  kt 

Table  2  Controlled  parameters  and  performance 
values  for  vertical  plane  altitude  change 


Similar  maneuvers  were  specified  for  flight 
tests  with  the  SR-71  [11]. 

After  each  simulation  run  the  test  pilot  an¬ 
swered  a  questionnaire  and  gave  a  Cooper- 
Harper  rating  [19]. 

4  Flight  Tests 

Dryden  Flight  Research  Center  operates  SR-71 
research  airplanes  (Fig.  9)  capable  of  sustained 
flight  at  Mach  3+  in  altitudes  above  80,000  feet. 
With  these  aircraft’s  flight  tests  for  similar  top¬ 
ics  as  considered  in  this  article  have  been  per¬ 
formed  [11]. 

In  the  high  supersonic  speed  regime  at  which 
the  SR-71  operates  flight  dynamics  similarities 


Fig.  9  Research  aircraft  SR-7 1  operated  by 
Dryden  Flight  Research  Center 


with  respect  to  hypersonics  exist.  For  example 
there  is  already  a  path  attitude  lag  of  4  to  5  sec¬ 
onds  at  Mach  3  so  that  an  alternate  vertical  speed 
feedback  control  technique  was  considered  [11]. 

5  Handling  Qualities  Results 

During  the  simulation  evaluations  the  short  term 
dynamics  were  kept  within  Level  I  [16]  limits. 

5.1  Path- Attitude  Lag  Variation 

Simulation  experiments  for  aircraft  configura¬ 
tions  with  different  path-attitude  lags  were  per¬ 
formed  using  the  “aggressive  stick”  piloting 
technique  for  controlling  flight  path. 

Handling  qualities  ratings  for  different 
path-attitude  delays  and  two  maneuvers  are  pre¬ 
sented  in  Fig.  10.  Despite  the  large  variation  of 
path  attitude  lag  there  are  only  small  changes. 
Vertical  plane  altitude  change  maneuvers  are 
rated  better  than  the  horizontal  turn. 

Similar  results  have  been  obtained  at  flight 
test  with  the  SR-7 1  research  airplane  at  high  su¬ 
personic  speeds  (Fig.  10). 

5.2  Effect  of  Vertical  Speed  Information 

The  “aggressive  stick”  technique  as  described 

above  requires  a  precise  indication  of  h  with  a 
high  resolution  and  in  an  adequate  format.  This 
was  done  with  the  diamond  symbol  moving 
along  the  pitch  ladder. 


Path-Attitude  Lag  [s] 

Fig.  10  Effect  of  path-attitude  lag  on  pilot  rat¬ 
ings  for  different  maneuvers  (GHAME 
model,  M=  10,  h= 1 10,000  ft) 
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Vertical  Speed  Indication  Time  Delay  [s] 
Fig.  11  Effect  of  time  delay  in  the  vertical 
speed  indication  on  pilot  ratings 
(GHAME  model,  M=10,  h= 1 10,000  ft) 


Time  Delay  of  Vertical  Speed  Information 
Results  concerning  the  effect  of  delayed  indica¬ 
tion  of  the  vertical  speed  (diamond  and  linear 
vertical  speed  and  acceleration  scales)  are 
shown  in  Fig.  1 1 .  There  is  a  significant  decrease 
in  pilot  rating  with  an  increase  in  time  delay.  At 
higher  values  the  pilot  was  only  barely  able  to 


use  the  “aggressive  stick”  technique.  Similar 
results  have  also  been  obtained  with  a  SR-71 
when  the  effect  of  an  altimeter  time  lag  has  been 
evaluated  (Fig.  11). 

Scaling  of  the  Diamond 

Fig.  12  presents  the  effect  of  the  diamond  scale 
on  handling  qualities  ratings  for  the  basic  simu¬ 
lation  model  at  Mach  10  in  110,000  ft.  In  this 
case  all  eigenvalues  are  stable  even  tough  the 
phugoid  has  a  time  to  half  amplitude  of  46,740 
seconds.  There  is  no  significant  effect  on  han¬ 
dling  qualities  ratings  (with  a  difference  of  1 
corresponding  to  normal  scatter). 

In  the  presence  of  a  phugoid  instability, 
there  are  significant  changes  (Fig.  13).  For  dia¬ 
mond  scales  with  a  resolution  of  less  than 
0.3  degrees  pitch  ladder  per  feet/second  vertical 
speed  the  handling  qualities  rating  gets  worse 
and  the  pilot  looses  control.  At  values  above  this 
threshold,  i.e.  at  a  higher  resolution,  the  rating 
remains  constant.  This  holds  for  two  different 
magnitudes  of  instability. 

With  respect  to  height  mode  instability 
there  is  not  such  an  effect  (Fig.  14).  Only  at 
diamond  scales  with  higher  sensitivity  there  is 
an  deterioration  of  the  handling  qualities  ratings. 


,  „  ,  r  deg  Pitch  Ladder  1 

D, amend  Scale  [Ws  Vert|ca|  Speed] 


Fig.  12  Effect  of  different  diamond  scales  on 
pilot  ratings  for  basic  GHAME  model 
(M=10,  h=\  10,000  ft) 


Diamond  Scale 


[deg  Pitch  Ladder  1 
ft/s  Vertical  Speed] 


Fig.  13  Effect  of  different  diamond  scales  on 
pilot  ratings  for  phugoid  instabilities 
(GHAME  model,  M=10,  h=l  10,000  ft) 
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Diamond  Scale 


[deg  Pitch  Ladder  ] 
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Fig.  14  Effect  of  different  diamond  scales  on 
pilot  ratings  for  height  mode  instabilities 
(GHAME  model,  M=10,  //=  11  (),()()()  ft) 


6  Conclusions 

In  high  speed  flight,  there  is  a  large  time  lag 
between  pitch  attitude  and  flight  path  responses. 
It  is  shown  that  there  is  an  upper  limit  for  this 
path-attitude  lag  at  a  high  level.  This  is  due  to 
the  coupling  of  two  real  zeros  in  the  pitch  atti¬ 
tude  transfer  function  to  complex  zeros. 

A  piloting  technique  (“aggressive  stick”)  to 
overcome  the  flight  path  control  problem  related 
to  the  path-attitude  lag  by  precise  flight  path  an¬ 
gle  or  vertical  speed  tracking  is  described.  This 
control  technique  also  has  an  stabilizing  effect 
on  the  phugoid. 

Simulation  experiments  show  that  there  is 
no  effect  of  a  path-attitude  lag  on  handling 
qualities  ratings  when  the  “aggressive  stick” 
technique  is  used.  To  utilize  this  piloting  tech¬ 
nique  a  precision  vertical  speed  indication  with 
appropriate  scaling  and  without  time  lag  is  nec¬ 
essary.  The  simulator  results  are  supported  by 
similar  results  obtained  from  flight  tests  with  a 
SR-7 1  research  aircraft. 
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